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Abstract 

We investigate the effect of exchange and correlation (xc) on the plasmon spectrum and the 
Coulomb drag between spatially separated low-density two-dimensional electron layers. We adopt 
a new approach, which employs dynamic xc kernels in the calculation of the bi-layer plasmon spectra 
and of the plasmon-mediated drag, and static many-body local field factors in the calculation of the 
particle-hole contribution to the drag. The spectrum of bi-layer plasmons and the drag resistivity 
are calculated in a broad range of temperatures taking into account both intra- and inter-layer 
correlation effects. We observe that both plasmon modes are strongly affected by xc corrections. 
After the inclusion of the complex dynamic xc kernels, a decrease of the electron density induces 
shifts of the plasmon branches in opposite directions. And this is in stark contrast to the tendency 
obtained within the RPA that both optical and acoustical plasmons move away from the boundary 
of the particle-hole continuum with a decrease in the electron density. We find that the introduction 
of xc corrections results in a significant enhancement of the transresistivity and qualitative changes 
in its temperature dependence. In particular, the large high-temperature plasmon peak that is 
present in the random phase approximation is found to disappear when the xc corrections are 
included. Our numerical results at low temperatures are in good agreement with the results of 
recent experiments by M. Kellogg et al, Solid State Commun. 123, 515 (2002). 

PACS numbers: 71.45. Gm, 73.20.Mf, 73.63.Hs 
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I. INTRODUCTION 



In a double layer two-dimensional electron system (2DES) the additional degree of free- 
dom associated with the layer index plays a key role in many recently discovered collective 
phenomena 1 . Suffice it to mention the freshly provided evidence that in the quantum Hall 
regime, near the half filling of individual layers, strong electron-hole inter-layer correla- 
tions cause Bose-Einstein condensation of excitons to occur in bi-layer 2D electron and hole 
systems^. In these structures, the ability to make separate electrical contacts to each layer 
has allowed the experimental detection of transport of neutral excitons as a counter-flowing 
current of electrons and holes in two layers. The independent control of layers is also critical 
for the experimental realization of the frictional drag^ which manifests itself under the 
condition that the second layer is an open circuit. A current driven along the first layer 
induces, via momentum transfer, a voltage in the second layer, which is measured experi- 
mentally. The drag effect complements and enriches the traditional experimental methods 
and in the last decade has proved to be a powerful tool to probe inter-layer electron-electron 
(e-e) interaction 5 . 

Recently the drag measurements have been extended to the limit of very low carrier 
concentrations. The dimensionless parameter r s = y/2/(kpa* B ), which describes the carrier 
density, n, and measures the strength of the electron-electron interaction^, varies approx- 
imately from 10 to 20 in the experiment on hole samples by Pillarisetty et alZ&. Here 
kp = v27m is the Fermi wave vector, and a* B = h 2 Ko/m*e 2 is the effective Bohr radius, with 
k,q the static dielectric constant and m* the effective mass of the carriers; a* B = 9.79 and 
3.45 nm, respectively, in the conduction (m* = 0.067mo) and valence {m* h = 0.19mo) bands 
of GaAs. In the experiment on electron samples by Kellogg et al. 9 , r s varies approximately 
from 2.5 to 4.5 and is appreciably smaller, owing to the small effective mass of the electrons. 
From a dimensional analysis it is clear, however, that in both types of samples, the Coulomb 
potential energy dominates the kinetic energy, and an adequate description of the drag can- 
not be provided by simple theories, which do not include the strong exchange-correlation 
(xc) effects. Consistent with this expectation, both experimental group o^-ff have found that 
the measured drag rate far exceeds (by more than a factor 200 in the hole samples) the value 
obtained from simple Boltzmann-equation calculations 10 . In addition, systematic deviations 
are observed from the quadratic temperature dependence of the drag rate (the characteristic 
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Fermi-liquid behavior of the Coulomb drag at low temperatures). 

One might hope that the above features could be understood in the framework of the 
phonon-mediated drag mechanis m n i 12 i 13 i 14 i 15 i 16 i 17 . However, the phonon-mediated drag pro- 
vides a dominant mechanism only in samples with large inter-layer spacing where the ab- 
solute value of the measured drag rate is much smaller than that observed in these experi- 
ments. Moreover, the phonon mechanism is not supported by the measured density-ratio- 
dependence of drag, which shows no peak at matched carrier densities in two layers, typical 
for that processes. Thus, both experimental group a 7 ' 8 ' 9 conclude that in these dilute 2DES 
the strong Coulomb interaction effects are mainly responsible for the observed new features. 

The carrier interaction effects on drag have been addressed previously in several exper- 
imenta l 18 i 19 i 20 and theoretica l 21 ' 22 ' 23 papers. While the theoretical prediction by Flensberg 
and Hu^ 2 - of a strong enhancement of drag by plasmons has been experimentally verified^ 9 - in 
high density electron samples, important differences have been reported 2 ^, from the results 
obtained within the random phase approximation (RPA). 

Lately, motivated by the recent set of low-density experiments^* 9 -, the drag resistivity has 
been calculated by Hwang et aZ.— and Yurtserver et ali^. Both works have included only the 
exchange interaction effects in the static limit via g-dependent but cu-independent local field 
factors (LFF). Besides, in the adopted approximation the nondiagonal inter-layer LFF have 
been taken to be zero^i, while for the intra-layer LFF the simple Hubbard approximation 
has been used, which significantly underestimates the LFF. As in the low density regime 
the inter-particle spacing in each layer becomes comparable with the inter-layer separation, 
one should expect that inter-layer correlations play an important role. At relatively high 
temperatures the dynamic xc effects are also critical. They can be especially important in 
the hole samples where, even at sufficiently low temperatures, the samples are effectively in 
the high temperature regime because of the small Fermi energy. 

In this paper we investigate the xc effects on the Coulomb drag in low-density e-e bi- 
layers. The spectrum of plasmons of the coupled 2DES and the temperature dependence of 
the drag resistivity are calculated in a broad range of temperatures. Our calculations include 
both intra- and inter-layer xc effects. We propose a new approach, which employs dynamic 
xc kernels in the calculation of the bi-layer plasmon spectra and of the plasmon-mediated 
drag, while the particle-hole contribution to the drag is still calculated by means of static 
many-body LFF. We find that the introduction of xc corrections results in a significant 
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enhancement of the drag resistivity and qualitative changes in its temperature dependence. 
In particular, a large high-temperature plasmon peak that is present in the RPA disappears 
when the xc corrections are included. At low temperatures our numerical results are in good 
agreement with the experimental findings by Kellogg et alA 

The paper is organized as follows. In Section II we outline the theoretical formulation 
of the frictional drag and provide the main formulas for its calculation beyond the RPA. 
In Section III a method is developed to relate the spin-averaged dynamic xc kernels and 
static many-body LFF of a bi-layer 2DES to the spin-resolved xc kernels and LFF of a 
single electron gas layer of the same total density. The results of our calculations are given 
in Sections IV- VI. First we present the spectra of bi-layer plasmons for different densities 
and discuss the effect of dynamic xc on the dispersion of each plasmon mode, comparing 
the dispersion curves with the corresponding curves obtained within the RPA. In Section 
V we study the temperature dependence of the drag resistivity at low densities within the 
RPA. Both plasmon and particle-hole contributions to the drag are considered. In the next 
Section we investigate the many-body xc effects on the drag resistivity for different densities 
and examine the distinctive features of the drag temperature dependence in comparison 
with the RPA. In the same Section VI we compare the obtained theoretical results with the 
experimental findings in low density e-e bi-layers. Finally, Section VII summarizes the main 
results of the paper. 



II. THEORETICAL FORMULATION 



We calculate the temperature dependence of the drag resistivity starting from the formula, 
obtained in the Kubo formalism^ 



h 2 dai do2 

Pd 



2e 2 a±a2TA dn\ dn^ 
Imn°(g,cj)Imn°(g,cj) 



(i) 



sinh 2 (^/2T) 

where A is the normalization area, hw and hlf the transferred energy and in-plane mo- 
mentum from the layer 1 to the layer 2, Wi2(q,u>) the dynamically screened inter- layer e-e 
interaction including the intra- and inter-layer many-body xc effects of a bi-layer 2DES, 
2 {q, oj) the finite temperature electron polarization function of an individual 2DES in the 
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absence of inter-particle Coulomb interaction. The imaginary and real parts of n? 2(^5 ^) 
are obtained from the Maldague formul a 22 ! 27 and are given by 
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Ren°(g,^) = g jl - e' 1 ' 1 - i- [M (t, £ (x, y)) 

- sign (C_ (a;, y)) M (t,C-(^3/))]} , (3) 

where we have introduced the following dimensionless quantities: t = T/ep, x = q/2kp, 
y = uj/Asfi n{t) = tin (e ly/ * — l), (±(x,y) = y/x±x with Bp and k F being the Fermi energy 
and the Fermi wave vector. The function 

POO 

F-i/ 2 (u) = l/^ dz/^(exp(z-u) + l) (4) 
Jo 

is the Fermi integral, the function 

M(t,u) = l/4t dzVu^z/ cosh [(z - fi(t)) /2t} 2 , (5) 
Jo 

and go = m*/irh 2 is the density of states of the noninteracting 2DES. Further we assume 
that in Eq. ([TJ the layer conductivities o\p, depend linearly on the electron densities n\^- In 
general, this is an acceptable approximation for the conducting 2DES in the magnetic-field- 
free case. Deviations from the linear regime can occur at very low electron densities. Here 
we consider only the regime corresponding to r s < 5; in this regime we expect the neglect 
of the corrections to the (7^2 oc n 12 to be a reasonable approximation. 

Our starting equation ((TJ is an approximation, strictly justifiable only to second order 
in the Coulomb interaction. The correct expression for the transresistivity involves a force- 
force response function, which is a four-point response function and cannot be reduced, in 
general, to a product of two-point response functions. In order to include the xc effects on the 
dynamically screened interlayer Coulomb interaction, in our calculations we have exploited 
an approximation scheme for W\2 proposed by Zheng and MacDonald^ (see also Ref. 23). 
Within this scheme an exact matrix Dyson equation for the Coulomb propagator—, W(q, 00), 
is approximated by 

W(q, u) = V e s(q, w) - V cS (q, u)IL(q, uj)V eS (q, u) (6) 



with the full polarization function, fl(q, u), defined in terms of the non-interacting polariza- 
tion function, fl (q,u), and unscreened effective Coulomb interactions, V e s(q,uj), as 

= — JM^ . (7) 

1 + V eS (q,UJ)Il (q,u) 

Here all quantities are 2x2 matrices. The effective Coulomb interactions V e s (whose explicit 
form is given below) build in Eqs. El and [7| the vertex corrections as sums of ladder diagrams 
evaluated within a "local approximation" 30 . Thus, the non-diagonal interaction matrix 
element Wi 2 (q, to) is given by 

W 12 {q,u) = — 8 

with the bi-layer screening function 

e hi (q, w) = £i{q, uj)e 2 {q, u) - V cS , 12 {q, w) 2 II?(g, uj)U° 2 (q, u) . (9) 

Here we introduce the screening functions of individual layers 

e 1>2 (q, u) = l + V cS>Ut22 (q, u)U l 2 {q, u) . (10) 

The approximation given by Eqs. (P)-()10j) accounts for the effect of dynamical screening 
of the Coulomb potential via the bi-layer screening function, eu(q, and includes the 
correlations between the test electron and the induced charge via the effective Coulomb 
interactions, V e g(q, uj). In this approximation xc LFF manifest themselves only via the 
intra- and inter-layer unscreened effective interactions 

Kff,y(<7, w) = v(q) (1 - G XCiij (q, u)) F^qd, qA) (11) 

where the intra- and inter- layer LFF G XC} ij(q,u)) decrease effectively the bare Coulomb in- 
teraction, v(q) = 2ire 2 /K q, by a factor of 1 — G XCj ij(q,uj) = 1,2 are the layer indices). 
Notice that the RPA is recovered if all the LFF in Eq. |H]are set to zero. 

Further we assume that the electrons in each layer are confined in a square well potential 
of infinite height and that scattering processes take place only within the lowest subband of 
each well. Assuming p±(z) = (2/cf) sm(nz/d) 2 and p 2 (z) = pi(z + A) as explicit forms for the 
electron density profiles in layers 1 and 2, respectively, one obtains^ for the corresponding 
intra- and inter-layer form factors (see Eq. fTTj) 

F , n _ 8^ 2 + ^ _ 32vr 4 (l - e-Q 



and 



„ , ^ 647r 4 sinh(ry/2) 2 , 
r\ A (Air z + rf) 



where A is the center-to-center inter-layer separation and d the width of quantum wells. 



III. DYNAMIC XC KERNEL AND STATIC LFF 

In general the local field factors of a bi-layer system, G^(q,u), depend on four indices: 
the spin indices \i and v, taking values j and j, and the layer indices i and j taking values 
1 and 2. But the LFFs that appear in Eq. ([lip, and hence in the effective interaction W% 2 , 
are averaged over the spin indices, at least as long as we assume that our bi-layer is not 
spin-polarized. Thus we have 

<Wff,a;) = \j2G%(q,u>) . (14) 

Unfortunately, reliable calculations of the LFFs for bi-layers (for example by the diffusion 
Monte Carlo method) have not been done to date. In order to proceed in the absence of this 
crucial information we propose an alternative approach, in which the spin-averaged LFFs 
of the bi-layer are expressed in terms of the spin-resolved LFFs of a single electron layer of 
the same total density n = n\ + n 2 - the "reference monolayer". We suggest to exploit the 
analogy between the spin index of the reference monolayer and the pseudospin (layer) index 
of the original bi-layer, neglecting the spin index in the bilayer altogether. We thus map the 
layer index (in the bilayer) onto the spin index in the monolayer and represent the intra- and 
inter-layer LFFs of the bi-layer in terms of the "spin-channel" and "charge-channel" LFFs 
of the reference monolayer, denoted respectively by G~(g, u) and G^(q,u) in the following 
manner: 

G XCill (q,u) = G+(q,uj) + G~(q,uj) , 

G xc , l2 (q,uj) = G+(q,uj) -G~(q,u) . (15) 

Of course this simple representation is possible only in the symmetric case, ri\ = n 2 , when 
we have G XC)U (g,u;) = G XC:22 (q, v) and G xc>12 {q,uj) = G XCj2 i(g, u). 

It is evident that Eq. ()15|) is not exact: for one thing, the interaction between electrons 
in different layers is different from the interaction between two electron in the same layer - a 
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difference that does not exist between parallel and antiparallel-spin electrons in the reference 
monolayer. Subtle differences persist even in the limit of zero inter-layer separation. In 
the reference monolayer the Pauli exclusion principle keeps two electrons of the same spin 
orientation from coming close to each other; but in the bi-layer two electrons in the same 
layer (i.e. with parallel isospin) can come together provided their real spins are opposite. 
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FIG. 1: Bi-layer intra- and inter-layer pair-correlation functions gi\(r) and §12(1*) and spin resolved 
pair-correlation functions g™(r) and g™(r) for a monolayer as directly obtained from the quantum 



Monte Carlo simulations, respectively, from Refs. 
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In spite of these difficulties we find that the approximation (|15|) is quite reasonable at 
small inter-layer separations. We get a good idea of the electronic correlations in a sys- 
tem by examining its pair-distribution functions. In Fig. Q we compare the spin-averaged 
pair-distribution functions of a bi-layer, gi\{r) (intra-layer) and g\ 2 { r ) (inter-layer), as ob- 
tained from the diffusion Monte Carlo simulations 31 , with the spin-resolved pair-distribution 
functions g^{r) and gV[(r) of the reference monolayer,— 2 . It is seen that there is an over- 
all good agreement between the corresponding distribution functions, i.e. g\\ agrees with 



gVl and g^ 2 with g*\. This agreement provides much support for the use of Eq. E3 As 
A increases, the agreement between the bi-layer pair-distribution function gfy (r) and the 
monolayer pair-distribution function g™(r) will deteriorate. However, the disagreement be- 
tween the interlayer correlations and the unlike spin correlation in the reference monolayer 
is more severe at small r than at large r^i. The pair-distribution function at small values of 
r, in turn, should control the behavior of the LFF at large wave vector q ^> 2kp, while the 
main contribution to the Coulomb drag, at low densities, come from scattering processes 
with q near 2kp. Thus, we see that our approach could be justified for the problem at hand, 
even when the inter-layer spacing is less than or comparable to the in-layer inter-particle 
distance, A < r s a* B , which is the case in the experiments of interest here. 

Let us now come to the discussion of the LFFs in the reference monolayer. Here we 
have to introduce further approximations. When treating the plasmon contribution to the 
Coulomb drag in our calculations we approximate the dynamic xc kernels, f XCti j(q,uj) = 
—v(q)G XCt ij(q,u) by their long- wavelength limits and employ the dynamic spin-spin and 
charge-charge xc kernels, evaluated by Qian and Vignale in Refs. and 0. In this way 
we include the frequency dependence of the dynamic charge-charge LFF, while their wave 
vector dependence is implicitly assumed to be linear so that the charge-charge xc kernel is 
independent of q. This is a reasonable approximation in view of the nearly linear behavior of 
the exact static LFF versus q in the single layer—. On the other hand, the spin-spin dynamic 
xc kernel is strongly dependent on q and diverges as q~ 2 at small wave vectors 36 . Thus, the 
intra- and inter-layer dynamic xc kernels in the long wavelength limit are represented as 

■A(u>) 



f XC;ij (q,co) = B + (to) + (25y-l) 



q 2 



+ B-(u) 



(16) 



where A(u) and B ± (u) are finite functions of ui and have been calculated to the leading order 
in the Coulomb interaction in Refs. QandQ. The authors of this paper have evaluated the 
four-vertex diagrams in the expansion of the response function of individual 2DES in the 
low-frequency/high-density limit. The interpolation formulas for Im A(u) and Im B ± (uj) 
have been constructed by combining the results for the imaginary parts with sum rules and 
previously known exact behaviors in the high-frequency limits. Thus, having the imaginary 
parts of A(uj) and 5 ± (a;), we obtain the corresponding real parts from the Kramers- Kroning 
relations. 

Figs. 121 and H3 show plots of the imaginary and real parts of the intra- and inter-layer xc 
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FIG. 2: Imaginary part (figure on the left) and real part (figure on the right) of the intra-layer 
dynamic xc kernel, f xc ^i(q,uj), of a bi-layer 2DES as a function of u>. The kernels are shown for 
the total r s = 2, 3, and 5 and for two different values of the wave vector q in the units of g^ 1 where 
go = m* /irti 2 is the density of states for the noninteracting 2DES. 
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FIG. 3: Imaginary part (figure on the left) and real part (figure on the right) of the inter-layer 
dynamic xc kernel, f XCj i2(q,uj), of a bi-layer 2DES as a function of uj. The kernels are shown for 
the total r s = 2, 3, and 5 and for two different values of the wave vector q in the units of g^ 1 
where go = m* /irh 2 is the density of states of the noninteracting 2DES. 

kernels, f xc ,ii(<l, and /xc,i2(<L oo), which we actually use in our calculations of the spectrum 
of bi-layer plasmons and of the Coulomb drag mediated by these plasmons. Notice that the 
kernels are calculated for the shown values of r s = 2, 3 and 5, which are obtained for the 
total number of electrons, n, in a bi-layer. The corresponding in-layer values of r i p~ layer , 
calculated for the n/2 number of electrons in individual layers, are r i p~ layer = 2y/2, 3\/2 
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and 5\/2. At finite frequencies, the term in Eq. EE related to A(u), dominates both in the 
short- and long- wavelength limits. The terms related to B^iu) make essential contributions 
to f xci j(q,u) only in the small uj region, where A(u) vanishes as u> 2 . 

The frequency- dependence of the LFF in the particle-hole continuum (PHC) region (finite 
wave vector, low frequency) is still largely unknown. For this reason, in evaluating the 
particle-hole contribution to the drag resistivity, we use the static limit of the LFF. We 
take advantage of the analytical expressions recently obtained for G~ c (q) and G+ C (q), LFF2&, 
which accurately reproduce the diffusion Monte Carlo data 3 - 3 -, as well as the exact asymptotic 
behaviors in the large and small q limits. Notice that both the dynamic xc kernel o 34 ' 35 
and the static LFF 3 —, which we have used in our calculations, are evaluated at T = 
for an ideal 2DES of zero width. At present there are no calculations of the LFFs of 
a two-dimensional system that take into account the finite width of the well and/or the 
temperature dependence. We expect, however, that the temperature dependence of the 
LFFs should play a relatively minor role in comparison to the temperature dependence that 
we explicitly include in the Fermi factors and the non-interacting polarization functions. 
Similarly, the effect of the finite width of the well should largely be taken care of by our use 
of the form factors Fy in Eq. El In Fig. 0] we plot the intra- and inter-layer LFF G XC) ij(q) 
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FIG. 4: Static many-body LFF of a double layer 2DES. The intra-layer, G X c,n (figure on the left), 
and the inter-layer, G XCj i2 (figure on the right), are plotted as a function of wave vector q for equal 
electron densities in the two layers corresponding to the total r s = 1, 2, and 3. The dotted curve 
shows the LFF Gh{x) in the Hubbard approximation. 

which we actually use in our calculations of the drag resistivity in the PHC region. For all 
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the values of r s shown in the figure, the intra-layer LFF G xc ^i(q) become significantly larger 
than the LFF in the Hubbard approximation, Gn(x) = x/\Jl + 4x 2 when x is close to unity. 
This makes the effective intra-layer interaction weak in the region near x — 1 and results 
in a weak in-plane screening effect. This effect becomes especially important in the regime 
of low density bi-layers where the main contribution to drag comes from the large-angle 
scattering processes with x ~ 1. In Fig.|3]we also see that although G^^iq) is smaller than 
Gx Ci n(g), the two LFF are of the same order of magnitude. 

IV. BI-LAYER PLASMON DISPERSIONS: EFFECT OF DYNAMIC XC 

Below we consider the symmetric bi-layer 2DES (ni = 712) with inter-layer spacing A = 28 
nm and quantum wells width d — 18 nm, which correspond to the experimental situation of 
Ref.S 

In this section we calculate the spectrum of bi-layer plasmons for the total r s = 2 and 
3 l^ r in-iayer _ 2^/2 and both within the RPA and beyond it taking into account the 

dynamic xc effects through the intra- and inter-layer xc kernels. In the symmetric systems 
the bi-layer screening function is represented as 

£ hi (q, u) = e+(q, w)e_(g, w) , (17) 

where 

e ± (q,u) = l + V ± (q,u)n (q,uj) , (18) 

with V±(q,u)) = V r e ff i n(g, u>) ± Kff,i2(<?, oj). The spectrum of the collective excitations is 
obtained from the zeros of the real part of the bi-layer screening function so the dispersion 
is obtained from the solution of the equation 

Re £±(q, i0 O p . a.c 

(?)) = • (19) 

The imaginary part of the screening function describes the damping of the collective modes. 
In this approach the imaginary parts of the dynamic xc kernels are finite at finite frequencies 
and we obtain that even at zero temperature the collective excitations are Landau-damped 
outside the T = PHC region. Also, we use the finite temperature forms for the polarization 
function U°(q,ui) which itself has a twofold effect on the plasmon spectrum. The finite 
temperature contribution to the imaginary part of U°(q,uj) outside the T = PHC region 

13 



enhances the damping of the collective modes, while temperature-induced changes in the 
real part of U (q,u)) modify the domain of existence of the collective excitations. 




0.0 0.5 1.0 0.0 0.3 0.6 0.9 
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FIG. 5: Bi-layer plasmon spectra at zero temperature (figure on the left) and at finite temperature, 
T = 0.5Tp (figure on the right) for the total r s = 3. The hatched area shows the PHC region 
at T = 0. The dispersions curves are calculated within the RPA (dashed curves) and beyond the 
RPA, taking into account the full intra- and inter-layer dynamic xc corrections (solid curves). 

In Fig.Elwe plot the plasmon spectra at zero and finite temperatures for the total r s = 3, 
based on the above Eq. The spectra are calculated within the RPA and beyond it 

taking into account the full intra- and inter-layer xc kernels. The upper-lying "optical" 
branches in both approximations exhibit the expected square-root ~ yfq dispersion in the 
long wavelength limit, while the lower lying "acoustical" branches have the linear ~ q 
dispersion. The inclusion of the dynamic xc kernels forces both the optical and acoustical 
branches to enter the PHC at substantially smaller values of q. At finite T the dispersions 
of collective excitations show the same behavior in the long wavelength limit while at the 
upper edge of the wave vector interval in which the collective excitations exist, the group 
velocity of each type of plasmon, independently of the approximation used, becomes infinite. 

Fig. |Hl shows that, as the electron density decreases, both the optical and the acoustical 
plasmon dispersions, calculated within the RPA, move away from the boundary of the PHC. 
We find, however, that the many-body effects beyond the RPA result in a new tendency, 
namely, the optical and the acoustical modes repel each other when r s increases so that the 
acoustical branch becomes closer to the boundary of the PHC while the optical branch is 
still pushed away from it. As we see below all these changes, induced by the dynamic many- 
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FIG. 6: Bi-layer plasmon spectra at zero temperature. The hatched area shows the PHC region at 
T = 0. The dispersion curves are calculated within the RPA (figure on the left) and beyond the 
RPA by taking into account the effects of intra- and inter-layer dynamic xc (figure on the right). 
The solid and dashed curves correspond to the values of the total r s = 2 and 3, respectively. The 
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dotted curves reproduce the plasmon dispersions from Ref. 
high density, n = 1.5 • 10 11 cm -2 , bi-layer samples with the inter-layer separation kpA = 4. 



body effects, are directly reflected in the temperature dependence of the plasmon-mediated 
Coulomb drag. 

V. DRAG RESISTIVITY WITHIN THE RPA 
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FIG. 7: Plasmon-mediated drag resistivity as a function of scaled temperature within the RPA. 
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First we present the RPA based calculations of the drag due to the exchange of plasmons 
in low-density bi-layers. In Fig. [7| the transresistivity (scaled by the factor (T/Tp) 2 ) is 
plotted as a function of temperature for five different in-layer densities corresponding to 
r in-iayer _ 2^ 2^/2, 3, 3\/2, and 5. Although the magnitude of the drag resistivity increases 
with r s , it is seen that within the RPA the qualitative behavior of the plasmon- mediated drag 
does not undergo strong changes when the electron density decreases, and it resembles the 
behavior at high densities^. The behavior of the transresistivity versus temperature exhibits 
an upturn at low temperatures near 0.2Tp and a plasmon enhancement peak at higher 
temperatures. At low temperature the transresistivity decreases slightly with decreasing 
density. In contrast to this, at high temperatures the height of the transresistivity peak 
increases approximately linearly with r % ^- la v er i n the range 2 < r l ™- la v er < 5. As seen 
in Fig. [7| both the upturn temperature of the plasmon-mediated drag and the position of 
the plasmon enhancement peak increase with r ™- la v er . This is consistent with the above 
discussion of the plasmon spectra, since in the RPA the plasmon dispersion curves for both 
modes move away from the PHC when the density decreases. Notice, however, that the 
absolute value of the upturn temperature and the position of the plasmon-enhancement 
peak (not scaled by Tp) show opposite tendency and decrease with r l ^- la v er . This is due to 
the fact that the Fermi temperature, Tp, decreases quadratically with r l ^- la v er . 

In Fig.|H]we plot the total transresistivity, due to exchange of both plasmons and particle- 
hole excitations against the part of the transresistivity that is due only to exchange of 
particle-hole excitations within the RPA. One can see that for r s = 2, 2y/2, 3, 3y/2, and 5 
when the density decreases, the contribution to the drag, made by particle-hole excitations, 
increases faster than the plasmon contribution to drag. However, the plasmon contribution 
remains essential in all cases and results in the plasmon enhancement peaks, clearly seen in 
the behavior of the scaled transresistivity versus temperature. 

The dotted and dashed-dotted curves in Fig. |H1 reproduce, respectively, the calculations 
by Jauho and Smith of the particle-hole contribution to the transresistivity based on the 
static screening approximation^ and by Flensberg and Hu of the plasmon-enhanced drag 2 ^ 
in the dynamic RPA. These curves correspond to the situation of the first drag experiment^ 
on high density samples with n = 1.5 x 10 11 cm -2 and with a large center-to-center inter- 
layer spacing, A = 37.5 nm. Notice that the plasmon enhancement peak is most pronounced 
in Fig. |H1 However, this strong effect is mainly caused by the large inter-layer separation 
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FIG. 8: Coulomb drag resistivity as a function of scaled temperature within the RPA for the 
in-layer r s = 2, 2^/2, 3, 3\/2, and 5 shown, respectively, from bottom to top. The solid curves 
correspond to the total drag, mediated both by plasmon and by particle-hole excitations, while the 
dashed curves represent only the particle-hole contribution. The dashed-dotted and dotted curves 
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are for the situation corresponding to the first drag measurements of Ref. I39t they are rescaled 
by a factor of 20, and represent, respectively, the total drag and the drag due to exchange of 
particle-hole excitations only. The latter is calculated within the static screening approximation. 

which exceeds significantly the effective spacing of the bi-layers in the experiment in Ref. 
for which the solid lines are calculated. As far as the plasmon-mediated drag depends more 
weakly on inter-layer spacing than the particle-hole contribution to drag 2 -, the plasmon 
enhancement peak is more pronounced for the dotted curves. 



VI. DRAG RESISTIVITY BEYOND THE RPA 

1. The effect of dynamic xc on plasmon-mediated drag 

In this subsection we include the intra- and inter-layer dynamic xc within the scheme 
(P)- (fTTj) and study their effect on the plasmon-mediated drag by comparing the obtained 
results with those of RPA-based calculations. In Fig. the scaled transresistivity as a 
function of temperature is shown for three equal density bi-layer systems with the total 
r s — 1, 2 and 3 corresponding to the in-layer r ™~ layeT = v^2, 2\/2 and 3^/2. It is seen 
that the upturn temperature in the scaled transresistivity essentially decreases when the 
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xc effects are included. Despite slight differences, in both cases of r f- layer = V2 and 2^/2 
the plasmons begin to contribute heavily to the inter-layer e-e interaction at temperatures 
less than 0.1T F . At r f- layer = 3V2 the effect becomes even stronger. Whereas, according 
to the RPA calculation, the upturn of the plasmon-mediated drag should be about 0.2Tp. 
This difference is explained by the changes in the spectrum of bilayer plasmons induced 
by the dynamic xc effects and discussed in the previous section. At low temperature the 
plasmon-mediated drag is mainly determined by the acoustical plasmons, whose frequency 
is reduced by the dynamic xc: this makes the acoustic plasmons easier to excite thermally, 
and causes their contribution to the Coulomb drag to be larger than in RPA. 

When temperature increases, the average energy of plasmons that mediate the drag in- 
creases. Hence, the plasmon damping, caused by dynamic xc effects, becomes stronger with 
a consequent reduction in the drag transresistivity. Thus, at temperatures around 0.37V the 
plasmon-mediated drag obtained within the RPA exceeds the drag which takes into account 
the xc effects. 
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FIG. 9: Dynamic xc corrections to the plasmon-mediated drag for the in-layer r l s n ~ la y er = ^2, 2\/2, 
and 3\/2- The scaled transresistivity vs temperature is shown within the RPA (dashed curves) and 
beyond the RPA taking into account the intra- and inter-layer xc effects (solid curves). 

At still higher temperatures the plasmon-mediated drag shows a broad peak both within 
the RPA and beyond it. The position of this peak is determined mainly by the position 
of the optical plasmons insofar as at such high temperatures the acoustical plasmons are 
heavily damped and merged into the PHC. As seen from Figs. El the peak position is at 
higher temperature at the lower density, and it changes slightly depending on the chosen 
approximation. These observations are again consistent with the behavior of the plasmon 
spectra in Figs. 03 and H3 The main contribution to the plasmon-mediated drag is made 
by the plasmons with energies fko ~ T while the scaled transresistivity peaks occur at 
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temperatures T <Tp. Therefore, as seen in Figs. El and |Hl at such small values of hw/AeF < 
0.25, the positions of the optical plasmons, and hence the positions of the peak in the 
plasmon contribution to drag in Figs. remain approximately unaffected. Notice also that 
at even higher temperatures the differences in the magnitude of plasmon-mediated drag in 
the different approximations diminishes with increasing T. 

2. The combined effect of dynamic and static xc on the total drag 

In this subsection we investigate the combined effect of the dynamic and static many- 
body xc on the total drag, mediated both by plasmon and by particle-hole excitations. 
In Fig. El the drag resistivity is plotted vs temperature within and beyond the RPA for 
the total r s = 1 (r* n ~ 

-layer = ^2). It is seen that in both approximations the particle-hole 
contribution to the scaled drag resistivity first shows a slight dip followed by a peak at higher 
temperatures. In the RPA the plasmon contribution to drag enhances this peak and shifts 
it to even higher temperatures so that the total transresistivity shows a pronounced peak 
approximately at the position where the plasmon contribution has a peak. As seen from 
Fig. El the introduction of the static exchange-controlled LFF increases the peak height of 
the particle-hole contribution to the drag and shifts the peak towards lower temperatures. 
On the other hand, the plasmon-mediated drag is moderately suppressed by the dynamic 
xc corrections as discussed in the previous subsection (cf. Fig. EJ). Thus, the resulting peak 
in the graph of the total transresistivity vs temperature becomes smaller at relatively low 
temperatures, while at high temperatures the drag rate shows a monotonic decrease in T. 

The described quantitative and qualitative differences in the behavior of the scaled total 
transresistivity within and beyond the RPA becomes more pronounced at lower densities. As 
seen from Fig. ^2 for the in-layer r i p- la v er = 2 a/2 and 3^/2 the total transresistivity beyond 
the RPA as a function of temperature shows no peak at all and this is in stark contrast to 
the peaked behavior of the transresistivity within the RPA. The disappearance of the large 
high-temperature plasmon peak results from the strong increase of the drag resistivity at low 
temperatures: we ascribe this to the fact that the contribution to drag, made by large-angle 
inter-layer scattering processes, becomes dominant when the many-body xc corrections are 
included. 

As shown recently in Ref. 0, the large- angle scattering processes with q ~ play an 
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FIG. 10: The drag resistivity vs temperature for the in-layer r l s n ~ la y er = ^q, The solid and dashed 
curves represent the total transresistivity and the separate particle-hole contribution to the drag 
resistivity with the full intra- and inter-layer xc corrections. The dash-dotted and dotted curves 
are, respectively, the total transresistivity and the particle-hole contribution to the drag calculated 
within the RPA. 
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FIG. 11: The full xc corrections to the total transresistivity for the in-layer r Y l ~ layer = 2\/2 and 
3\/2- The dashed and solid curves represent the scaled transresistivity vs temperature, calculated, 
respectively, within and beyond the RPA. The RPA data are multiplied by a factor of 3 and 6, 

, ■ i p in— layer n /n j in— layer fK 

respectively for r s = 2v I and r s = 6\/2. 

important role in the Coulomb drag in low density bi-layers. In such samples with sufficiently 
small inter-layer separations, kpA is so small that the usual exponential cutoff, exp(-qA), of 
the inter-layer e-e interaction is no longer effective in suppressing the large-angle scattering 
processes. Thus, according to previous predictions^., the large-angle inter-layer scattering 
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events, due to the divergence of the scattering phase-space near q ~ 2k f, lead to a T 2 lnT 
behavior of the drag resistivity at low temperatures, in lieu of the usual T 2 temperature 
dependence, which is driven by small-angle scattering processes. To distinguish clearly the 
contributions to the drag made by the small- and large-angle inter-layer scattering processes, 
in Fig. ^1 we have plotted the dimensionless drag intensity, I(x), within and beyond the 
RPA versus the average transferred momentum, x = q/2kp, at low temperature T = 0.057V . 
We have normalized the drag intensity to its corresponding peak values, p , so that the 
transresistivity p D = p p f I(x)dx. It is evident that the xc effects at the low density reduce 
essentially the component from the small-angle scattering processes so that the large-angle 
scattering contribution dominates the drag. Thus, the drag rate is strongly enhanced at low 
temperatures owing to the logarithmic corrections provided by the large-angle scattering 
events. By contrast, in the RPA the contribution of the small-angle scattering processes 
remains dominant down to the lowest temperatures (see Fig. 112)1 . 




FIG. 12: Drag intensity I(x) as a function of transferred momentum x = q/2kF for a low density 
n = 3.8- 10 10 cm -2 bi-layer. For comparison the drag intensity is plotted also for a high density n = 
1.5 • 10 11 cm -2 bi-layer. The dashed (dotted) and solid (dash-dotted) curves represent respectively 
the results of calculations within and beyond the RPA at T = 0.05Tp for n = 3.8 • 10 10 cm -2 
(n = 1.5 • 10 11 cm" 2 ). The respective values of the drag resistivity are given by the area under 
these curves. The curves are normalized to their particular peak values, p RPA = 0.57 (0.06) and 
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2.71 (0.01) Q for n = 3.8 • 10 10 cm" 2 (n = 1.5 • 10 11 cm 



Lack of experimental measurements of the Coulomb drag at high temperatures does not 
allow for the time being an experimental verification of our predictions on the position 
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and strength of the plasmon peak in low density bi-layers. Notice however that at low 
temperatures, as shown in Fig. El our numerical results are in good agreement with the 
experimental findings by Kellogg et a/.— for n = 3.8 x 10 10 cm -2 . In this figure we have 
shown the drag resistivity due to exchange of particle-hole excitations. Inclusion of the 
dynamical part will only slightly increase the transresistivity since at such low temperatures 
plasmons only start to contribute. As seen in Fig. E3 the RPA calculations near T = 5 K 
underestimates the experimental results more than three time. Meanwhile after inclusion 
of the xc corrections, the drag resistivity is about 88% of the experimental result. This is 
very good agreement against the background of the usual discrepancy between theory and 
experiment in determining the magnitude of the drag. 
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FIG. 13: The drag resistivity versus temperature for n = 3.8T0 10 cm -2 . The solid curve represents 



calculations when the intra- and inter-layer many-body corrections are inc 
is the RPA results. The symbols are the experimental findings from Ref. 
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In Fig.^lwe plot the drag resistivity versus temperature for six different densities. From 
the comparison with the respective Fig. 1 from the experiment in Ref. it is seen that 
our results reproduce very well all the experimental findings for the densities n = 5.2, 4.7 
and 3.8 • 10 10 cm" 2 . For lower densities n = 1.7, 2.3 and 3.1 • 10 10 our calculations start 
to deviate from the experimental results. This is clearly seen also in Fig. ED where we 
show a log-log plot of the drag resistivity versus density n for three different temperatures: 
T — 1, 2, and 4 K. The symbols are our calculations while straight solid line represents 



n 4 dependence. This figure corresponds to Fig. 2 from Ref. where it has been found that 



the experimental results can be very well approximated by the quartic density dependence. 
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As seen in Fig. [T31 except very low densities, our calculations of the drag resistivity, which 
include the xc corrections, are also approximated by the quartic density dependence quite 
well. This contrasts the RPA-based results where the drag resistivity shows n~ 3 dependence. 




2 4 

7[K] 

FIG. 14: The drag resistivity vs temperature for six different densities corresponding to Fig. 1 from 
Ref. 9. The densities are in units of 10 10 cm~ 2 . 

VII. SUMMARY 

In this paper we have investigated the role of the many-body xc effect on the dispersion of 
collective excitations and on the frictional Coulomb drag in low density bi-layer 2DES. We 
have calculated the spectrum of the bidayer plasmons and the temperature dependence of the 
drag resistivity in a broad range of temperatures. Our calculations include both intra- and 
inter-layer xc effects. These effects are critical in low density bi-layers since the inter-layer 
separation in these structures is comparable with the inter-particle spacing in an individual 
layer. We have proposed a new approach and employed the full dynamic xc kernels of the 
bi-layer 2DES in calculations of the acoustical and optical plasmon dispersions as well as in 
the drag, mediated by exchange of bi-layer plasmons. In studying the many-body effects on 
the particle-hole contribution to drag we have still used the static exchange-controlled LFF. 
The combined effect of the dynamic and static many-body xc on drag has been investigated 
in comparison with the RPA-based calculations of the transresistivity vs temperature for 
different values of the electron density. At low temperatures we have compared our results 
with the experimental findings from Ref. 0. 

We have observed that after the inclusion of the full dynamic xc kernels, a decrease of the 
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FIG. 15: Log-log plot of the transresistivity vs density for three different temperatures: T = 1,2, 
and 4 K. This figure corresponds to Fig. 2 from Ref. |9|. The solid lines are proportional to n~ 4 . 



electron density induces shifts of the plasmon branches in opposite directions. And this is 
in stark contrast to the tendency obtained within the RPA that both optical and acoustical 
plasmons move away from the boundary of the PHC with a decrease of the electron density. 
In the presence of the xc effects, the optical mode still moves away from the PHC while 
the acoustical branch moves closer to the PHC. The effect is so strong that at the total 
r s = 3 the acoustical branch enters the PHC and becomes completely destroyed at finite 
temperatures. These changes of the plasmon spectrum, induced by the intra- and inter-layer 
dynamic xc, are reflected in the temperature dependence of the plasmon-mediated drag. 

After inclusion of the xc effects the upturn temperature, which is determined by the ther- 
mal excitations of the acoustical modes, exhibits a "redshift" while the plasmon enhancement 
peak, which is mainly due to the thermal excitation of high energy optical modes, under- 
goes to a "blueshift" with decreasing carrier concentration. Already for the total r s = 2 the 
upturn temperature decreases so strongly that the bilayer plasmons heavily contributes to 
the drag starting at temperatures less than O.lTp, which is approximately a factor of 2 small 
than the upturn temperature calculated within the RPA. The effect becomes even stronger 
for the total r s = 3 so that the upturn temperature disappears practically. Our results show 
that at low temperatures the dynamic xc effects enhance the plasmon-mediated drag while 
at high temperatures they have the opposite effect. We ascribe this to the enhancement of 
the damping of high energy plasmons, caused by dynamic many-body xc effects. 

The combined effect of the dynamic and the static many-body xc in low density bi-layers 
results in both quantitative and qualitative changes in the behavior of the drag resistivity. 
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At temperatures near 0.5Tp the transresistivity with the xc effects included is larger approx- 
imately by a factor of 3 than that obtained within the RPA for the in-layer r l ^- la v eT = 2y / 2. 
The enhancement increases strongly at low temperatures where the drag is dominated by 
large-angle scattering processes with q ~ 2kp- We have found that the plasmon enhancement 
peak, which is clearly seen in the RPA based calculations, disappears at large r s after the 
full xc corrections have been introduced. The disappearance of the large high-temperature 
plasmon peak is the result of the strong enhancement of the particle-hole contribution to 
drag at low temperatures and of the moderate suppression of the plasmon-mediated drag 
at high temperatures when, respectively, the full static and dynamic xc corrections are in- 
cluded. We have ascribed this behavior to the change in the specific contributions to drag 
made by small- and large-angle scattering events. Our calculations demonstrate clearly 
that the xc effects increase the large-angle scattering component of the drag making it the 
dominant contribution. This component strongly enhances the drag for two reasons. First, 
at low temperatures the e-e scattering phase-space diverges near x = q/2kp — 1. On the 
other hand, when x ~ 1, the static intra-layer LFF Gn(x) becomes close to unity, lead- 
ing to a reduction of the effective intra-layer interaction V^nix). This by itself weakens 
the dynamic intra-layer screening and enhances the drag. Thus, against the background of 
the large transresistivity at small T, the plasmon-mediated contribution to drag does not 
result in a peaked behavior of the scaled transresistivity and the total drag rate decreases 
monotonically in T. In the RPA, the drag still contains a strong small-angle scattering 
component which suppresses partially the enhancement at small T and the plasmon peak 
remains visible. 

Currently, the available experimental data for the low density electron samples are re- 
stricted to low temperatures where the plasmon contribution to drag is small. An experi- 
mental verification of the plasmon peak disappearance requires new drag measurements in 
low density electron bi-layers in a broad range of temperatures close to Tp. 

At low temperatures T <C Tp our numerical results are in good agreement with the ex- 
perimental findings by Kellogg et al. 9 . At the upper edge of temperature interval where 
the measurements are available our calculations of the drag resistivity give about 90% of 
the measured drag. This agreement is especially important against the usual discrepancy 
between theory and experiment in determination of the drag absolute value. Our calcula- 
tions beyond the RPA show that the drag resistivity as a function of density exhibit n~ A 
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dependence for not very low densities. This is in contrast to the RPA based calculations of 
the drag resistivity and in agreement with the experimental findings. 
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